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Abstract. For non-abelian simple objects in a unitary modular category, the 
density of their braid group representations, the #P-hard evaluation of their 
associated link invariants, and the BQP-completeness of their anyonic quantum 
computing models are closely related. We systematically study such properties 
of the non-abelian simple objects in metaplectic modular categories, which are 
unitary modular categories with fusion rules of S0(m)2 for an odd integer m > 3. 
The simple objects with quantum dimensions ^Jm have finite image braid group 
representations, and their link invariants are classically efficient to evaluate. We 
also provide classically efficient simulation of their braid group representations. 
These simulations of the braid group representations can be regarded as qudit 
generalizations of the Knill-Gottesmann theorem for the qubit case. The sim- 
ple objects of dimension 2 give us a surprising result: while their braid group 
representations have finite images and are efficiently simulable classically after 
a generalized localization, their link invariants are #P-hard to evaluate exactly. 
We sharpen the #P-hardness by showing that any sufficiently accurate approx- 
imation of their associated link invariants is already #P-hard. 



1. Introduction 

Unitary modular tensor categories (UMCs) are intricate algebraic structures 
arose in a variety of fields, in particular in the study of topological quantum field 
theories (TQFTs) [Tuj . conformal field theories (CFTs), and topological phases 
of matter. Mathematical constructions of UMCs include representation theories 
of quantum groups at roots of unity and vertex operator algebras. In condensed 
matter physics, each UMC is a theoretical anyonic system, and a simple object 
models an anyon. Therefore, we will use the words simple object in a UMC and 
anyon interchangeably. In topological quantum computation, any anyon non- 
abelian can be used to construct an anyonic quantum computing model |Waj . 
The realization of UMCs in real physical systems and the universality of anyonic 
quantum computing models inspire many new mathematical problems. 
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2 METAPLECTIC UMCS 

In this context, the most studied sequence of UMCs is SU{2)k,k = 1,2,---. 
One reason is the conjectured relations between SU(2)2, SU(2)4, and SU(2)3 and 
the potential non-abelian statistics in fractional quantum Hall states with filling 
fractions u = |, |, and y |RR] . By the level-rank duality, this sequence is es- 
sentially SU(A;)2 — the A-series at level two. In this paper, we will focus on the 
5-series at level two — S0(m)2 for an odd integer m > 3. (The level-rank duality 
for S0{m)2 is more involved, in particular SO {171)2 is not the same as S0{2)m)- 
We will call any UMC with the same fusion rules of SO {171)2 for some odd integer 
m > 3 a metaplectic modular category. Our choice of metaplectic is motivated 
by the connection of such UMCs and the metaplectic representations of the finite 
symplectic groups Sp{2n, F.^) when m is a prime p. 

There are several reasons to be interested in metaplectic modular categories. 
One is their conjectured property F and related explicit locality, while the other is 
the solutions that they provided for the generalized Yang-Baxter equations |RW[ 
IKWj . But our interest mainly comes from their potential realization in condensed 
matter systems. With this application in mind, we study problems inspired by 
their application to quantum computing, in particular analyzing their potential for 
topological quantum computing. In the companion paper [HNWj . we addressed 
their relevance to some proposed physical systems. 

Let X be a simple object in a UMC C. We call X non-abelian if its quantum 
dimension dx > ^ because X will model a non-abelian anyon in physics; otherwise, 
X is abelian, and then dx = ^- We are mainly interested in non-abelian simple 
objects. Associated to each simple object X in a UMC C is a unitary braid group 
representation px and an isotopy invariant 'Jx{L) of framed links L [Tuj . Modulo 
subtleties of encoding, the representation matrices of px can be used as quantum 
circuits for an anyonic quantum computing model Ax |Wa] . We are interested 
in the BQP-completeness of such anyonic quantum computing models. Then the 
density of the braid group representations px, the #P-hardness of evaluating 
'Jx{L), and the BQP-completeness of Ax are all closely related. 

For the non-abelian simple objects X of quantum dimensions dx = 2cos{-^) 
in SU(2)fc, the density, the #P-hardness, and the BQP-completeness all match 
perfectly |FLWH rFLW2j . Non-abelian simple objects with all three properties will 
be called strong, and otherwise called weak. Note that abelian simple objects are 
always weak in this sense. The weakest non-abelian object will have finite image 
braid group representation, polynomial time algorithm for the evaluation of the 
associated link invariant, and a classically efficiently simulable computing model 
Ax- In the sequence SU(2)fc, non-abelian simple objects of quantum dimensions 
dx = 2cos(^^) are weak if and only if k = 2,4. Note that the case k = 1 
is abelian. When k = 2, this is the Ising anyon a, the finiteness of the braid 
group representation, polynomial time computable link invariant, and classically 
simulable computing model are all known |Jo83t IJo87j . So the Ising anyon is weak 
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in all three aspects, where the classical simulation of the model follows from 
the Knill-Gottesmann simulation of Clifford circuits (See e.g. |NCj ). For level 
k = 4, the finiteness of braid images and polynomial computable invariants are 
also known |Jo87] . In this paper, we show that efficient classical simulation of 
braidings is also possible after a localization. Therefore, a natural question is if it 
is possible for a non-abelian anyon to be weak in only one or two aspects. 

We systematically study non-abelian simple objects in metaplectic modular cate- 
gories. We prove that the simple objects with quantum dimensions ^/m in SO {711)2 
are weak in all three aspects, but the simple objects of quantum dimension 2 in 
JKr=6 give us a surprising result: while their braid group representations have fi- 
nite images and are efficiently classically simulable after a generalized localization, 
their link invariants are #P-hard to evaluate exactly. Actually we prove that any 
sufficiently accurate approximation of their associated link invariants is #P-hard. 

The contents of the paper are as follows. In Section 2, the basic data for meta- 
plectic modular categories are listed. In Section 3, we analyze the braid group 
representations associated to the quantum dimension 2 simple objects. It follows 
from our analysis that their braid group representations have finite images. In 
Section 4, we provide classical efficient simulation of the Gaussian braid repre- 
sentations and the braid representations from |KWj . These simulations can be 
regarded as qudit generalizations of the Knill-Gottesmann theorem. In Section 5, 
we prove that any sufficiently accurate approximation of the link invariant associ- 
ated with the quantum dimension 2 simple objects is already #P-hard. 

2. Metaplectic Modular Categories 

One systematic way to construct UMCs is via the representation theory of quan- 
tum groups at a particular root of unity |Tu] . For each simple Lie algebra g and 
an integer k > 1, called the level of the theory, two UMCs can be constructed 
from the representation theory of Ug{Q), where q = e^T^,/ = mg(k + kg), where 
hg is the dual Coxeter number of g and rrig = 1 for A, D, E, 2 for B, C, F4, and 3 
for G2- For SU(2), /isu(2) = 2 and for SO(m), hso{m) = m — 2. The two UMCs for 
g = e^T are conjugates of each other, so for definiteness Gk denotes the category 
from q = , where G is the simple Lie group with Lie algebra q. 

2.1. S0(m)2 UMCs. The UMC S0(m)2 is of rank r + 4, where m = 2r + 
l,r > L Their simple object representatives will be denoted as Irr{S0(m)2} = 
{1, Z, X^, X^i, Yj, 1 < j < r}. In |NRj . 1 is denoted as Xq, Z as X2A1 , and Yj as Xy . 
Their quantum dimensions are di = dz = I, dx^ = dx^, = dy^ = 2,1 < j < r. 

All fusion rules of those UMCs can be deduced from the ones listed below. In 
the following, we also denote 1 by Yq sometimes. 

(1) X,^X, = (B'jZlYj 

(2) X, ® Yj = X,®X'^,l<j <r 
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(3) x,®x'^ = z® ©^=;y,- 

(4) Z®X, = X[ 

(5) Z®Z = 1 

(6) Z®Yj = Yj, 1 < J < r 

(7) Yj®Yj = l®Z ® Fmin{2i,m-2j}, 1 < J < r 

(8) ® Yj = Y\i_j\ ® ymin{i+j-m-i-j}, 1 < j < ^ 7^ j- 

Note that the fusion rules for the sub-category consisting of {1, Z,Yj,l < j < r} 
are exactly the same as those of the representations of the dihedral group Dm of 
order 2m. 

The twist or topological spin of a simple object X is of the form 6x = e^'^''*-^, 
where the rational number hx, only defined modulo 1, is called the scaling dimen- 
sion of X. For 5*0 (m) 2, we have: 

hz = 1, hx, = -, hx'^ = hy^ = , 0<J <r. 

Some braidings for S0(m)2 are as follow: 



_ g7ri(m+l)/m ^^I'^l = ^iri/m J^^l'^i = g7ri(m-l)/m 

2.2. Metaplectic modular categories. 

Definition 2.2. Lei C 5e a UMC, and X be a simple object. 

(1) C metaplectic if its fusion rules are the same as those of S0{m)2 for 
some odd integer m > 3. 

(2) X is called non-abelian if its quantum dimension dx > 1; otherwise, X is 
called abelian, and then dx = ^ 

(3) An abelian simple object is called a boson if its twist Ox = ^, a fermion if 
Ox = —1, and a semion if Ox = ±i. 

We will use the same notation for simple objects in S0{m)2 as for all metaplectic 
modular categories. The fusion of an abelian simple object with other simple 
objects only permutates them. From the above fusion rules, we see that the simple 
object Z of any metaplectic modular category is a boson. A complete classification 
of all metaplectic modular categories is open. 

When m = 3, metaplectic modular categories are of rank=5. In |BNRW] . 
all rank=5 UMCs are classified. In particular, there are 4 metaplectic modular 
categories for m = 3: the SU{2)4 and its conjugate, the Jones-Kauffman theory 
at A = ie~^ , denoted as JKr=Q, and its conjugate. The SU {2)4 and the JKr=e 
are distinguished by the Frobenius-Schur indicator of X^: —1 for SU{2)4 and 1 for 

JKr=6. 



(2.1) 
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The full data of S'[/(2)4 and JKr=% can be found in |Bo] and |KL] . respec- 
tively. While the F-matrices in |Boj are unitary, the F-matrices in |KLj need to 
be renormalized as in |Wa] to be unitary. 

Each simple object in a metaplectic modular category has an integral squared 
quantum dimension. For braid group representations, the afforded representations 
by X'^ and differ only by a one-dimension representation, so we will only discuss 
the one from X^. Also due to the similarities among Yi's, we will focus only on Yi. 
The equivalence of the braid representation of X^ in SO {171)2 and the Gaussian 
representation below is proven by E. Rowell and H. Wenzl (Roj . 

3. Extra-special p-groups and metaplectic representations 

3.1. Braid group representations. Suppose C is a UMC, and Irr{C} = {Xi}i^i 
is a set of simple object representatives, i.e. one from each isomorphism class of 
simple objects. Then each simple object X in C gives rise to a unitary represen- 
tation px = {px,n} of the braid group 23 = Such representations for the 
n-strand braid group !B„ can be conveniently described using an algebra: An^x = 
Hom(X®", X®"). The algebra An^x is an analogue of a regular representation for a 
finite group, and An^x — ©ie/ End(Vx®n xj, where Vx®-^,Xi = Hom(X®", Xj). The 
representation Vx®n of denoted also as px.n, can be conveniently described 
using graphical calculus. An orthonormal basis of Vx®-^,Xi is given by admissible 
labelings of any connected uni-trivalent tree [Waj . 

In |RWt IGHRj . localization and generalized localization are studied for a simple 
object X whose squared quantum dimension is an integer. A (generalized) localiza- 
tion of the braid group representation px is a single unitary (generalized) i?-matrix 
Rx such that the natural braid group representation from Rx-, denoted as pR_y, 
is equivalent (as a representation) to a braid representation px constructed from 
px as follows: decomposing each representation px,n of the n-strand braid group 
into its irreducible summands: px,n = ®jPx,n,j, then px,n = ®j^x,n,jPx,n,j for 
some non-zero multiplicities mx,n,j- 

The localizing qudit representation pji^ has a natural tensor product structure, 
similar to the quantum circuit model. Therefore, we will refer to a localizing 
representation pji^ of a simple object X as its qudit representation associated to 
Rx- 

For each such object X localized by a (generalized) i?- matrix Rx, we have 
two braid group representations px and Prx, which have the same irreducible 
summands. The question that we are interested in is: 

For each n, is the image pxi'Bn), or equivalently p^^{'Bn), a finite group? If so, 
which finite group ? 

3.2. Finiteness for the representation px^- 
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3.2.1. Extra-special p- groups. Let Xj (or Zi) be the Pauh matrices ax (or az) act- 
ing on the i-th qubit of the n-qubits (C^)®". The real Pauh group CP„ on ra-qubits is 
the group generated by aU such Pauh gates. (The complex Pauli groups are slightly 
more complicated.) As an abstract group, it is an extra-special 2-group isomorphic 
to the central product of n-copies of the dihedral group -D4 of order 8. Its outer au- 
tomorphism group Out(y„) = Aut(CP„)/Inn(y„) is an orthogonal group 0^(2n, F2) 
over the field F2 [Grj . The definition of the orthogonal groups 0^{2n, ¥2) does not 
use the usual bilinear forms, which would lead to the symplectic group Sp(2n, F2). 
Instead they are defined using quadratic forms over the field F2. Over the F2- 
vector space F^"", there are two inequivalent non-degenerate quadratic forms with 
maximal isotropic subspaces of dimensions n and n — 1, respectively. The orthog- 
onal groups 0^(2n, F2) are the linear transformations of F^" that preserve these 
two quadratic forms, respectively. For the Pauli group CP„, the quotient CP„/Z(CP„) 
is Ff" with a natural quadratic form q: q{v) = f ^ G Z(7n) — ^2,v G Fg". 

The Clifford group C„ on n-qubits is the normalizers of the Pauli group IP„ in 
the unitary group U(2") up to sign. The two groups fit into the following short 
exact sequence: 

l^T„^e„^Sp(2n,F2)^l, 
where Sp(2n, F2) is the symplectic group of the F2-vector space Fg". 
For an odd integer m > 3, let E^^_^ be the group generated by Ui, . . . ,m„_i 
with relations: 

(3.1) uT = 1 

(3.2) MjMj+l = ZUi+iUi 

(3.3) UiUj = UjUi, 1^— j| > 1, 

where 2; is a central element of order m. 

When m is an odd prime p and n > 1 is odd, then -Ep„_i is an extra-special 
p-group of exponent p. In physics, -Ep ^.i is often cahed the finite Heisenberg 
group. Then the outer automorphism group Outi{Ep „^_^) of Ep „^_^ (the subscript 
/ of Out J means that we consider only automorphisms that are the identity on the 
center of Ep^^_^) is Sp(?2 — 1, Fp) [Wi]. We have another split short exact sequence: 

1 ^ E;„_i ^ ep,„_i ^ Sp(n - 1, Fp) -> 1. 

Therefore, the extra-special p-group E^ ,^_^ can be regarded as the p-generalization 

of the Pauli group CP™^, while the semi-direct product of the extra-special n-group 

2 

by the symplectic group Sp(?T, — l,Fp), denoted be regarded as the 

generalization of the Clifford group Cn-i . 

2 

A similar interpretation can be made for n even. In this case, the center of 
^p,n-i is ® rather than Zp (the extra Zp is generated by U1U3 ■ ■ ■ Un-i-) We 
call this group an almost extra-special p-group, though this terminology is not 



METAPLECTIC UMCS 



7 



standard. Below we will see that those groups naturally arise as images of braid 
group representations. 

3.2.2. Gaussian representations . Let w be a primitive m-th root of unity and con- 
sider the C-group algebras ES{u, n — 1) = C[E^ J, where z = u~'^. Notice that 
ES{u,n — 1) has dimension m"~^ and is semisimple. Therefore, ES{u,n — 1) = 
©jM„.(C), where Mjj. (C) is the full Ui x Ui matrix algebra and X^i^? = m"'~^. 

The Gaussian representation |Jo89j 7 : S„ — )■ ES{u, — 1) is defined on braid 
generators of S„ by 

^Jm ^ — ^ 

Direct computation shows that this leads to a unitary representation of the 
braid group. To match the representation px, of S0{m)2 exactly, we need to add 

a phase factor e ^"s ' to 7(crj). Since this phase factor is not important for our 
discussion below, we will omit it for convenience. 

There is another representation of 23„ in ES{uj,n — 1): 

m-l 

m ^-^ 

j=0 

where t + t^^ + 2 = m. This representation is usually called the Potts represen- 
tation. While for m = 3, the Potts and Gaussian representations coincide, they 
differ for m > 5. 

3.2.3. Metaplectic representations . It was shown in |GJoj that the images of the 
Gaussian braid group representations are finite groups. In fact, for n odd the 
analysis in |GJo] shows that, projectively, 7(S„) is isomorphic to the finite simple 
group PSp(n — l,Fp) when m is a prime p. More carefully, we have: 

Theorem 1. When m is a prime p, the images 7(2^) of the Gaussian represen- 
tations are: 

(1) when n is odd, then 7(S„)/Z(7(!B„)) = Sp(n — l,Fp), where the center 
Z(7(S„)) ^ Z4 ifp = 3 mod 4 and Z2 © Z2 if p = 1 mod 4. 

(2) when n is even, then 7(S„) fits into the following exact sequence: 

1 ^ ^ 7(S„) ^ Sp(n - 2, Fp) ^ 1. 

For the braid group representation associated to the simple objects of S0{m)2, 
let Hn,x^ denote the images pxc,n{'S>n) of the "Bn from the object X^. Then for m 
a prime p, -ff^.x^ — 7(Sn) for all p [Roj . 

A uniform way to understand the images 7(S„) is to treat the short exact 
sequence for n even as the definition of the symplectic group Sp(n — l,Fp). The 
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image of the braid group in the Gaussian representation can be understood as 
follows. We can verify that: 

[7(o"i-i)]^Mi7((Ti_i) = uu^^Ui 
[7(o-i)]"^Mi7(ai) = Ui 
(3.4) [-/{aj)]^ Ui-f{aj) = Ui , \i - j\ > 1 

Braiding transformations are, therefore, automorphisms of Ep^^/Z{Ep^^). Hence, 
the image of the braid group is a subgroup of the group of automorphisms of 
Ep^^/Z{Ep,^). This is equal to the metaplectic representation [GJo] of Sp{n — 1, ¥p). 

3.3. Finiteness for the qubit representation pi?^^. Let Ry^ be the following 
8x8 matrix, which is the block sum of two 4x4 matrices: 



/ucos(^) 


isin(^) 








where u 







COSI — 



ism^^j 


z/cos(— ) 




V 





cos(^) 

^ m ' 










isml — , 



COS - 

^ m ' 



- i sin(- 

^ ri 





-1 if m = 3 and u 



We will call 





cos(^ 


-isin(^)y 

1 if m > 5. 
For the braid generator cTj G S„, set 

then pRy^ is a representation of !B„ on the {n + l)-qubit (^c2^(g)(n+i) 
this representation the qubit representation. This representation is a generalized 
localization of the py^ for Yi of the JKr=e and its generalizations to m > 5 [RWj . 
We do not know if there are localizations of py^ 's. 

By the computational basis of n-qudits (C")®", we mean the basis consists 
of the tensor products {\ji) (g) \j2) ® • ■ ■ (8) |jn),0 < ji,j2, ■ ■ ■ ,jn < m - 1} of 
the basis {\j),j = 0,..,m — 1} of C". When m = 2, we have a one-to-one 
correspondence between the computational basis of n-qubits and the n-bit strings 
|xi ■ ■ ■ Xn),Xi G {0, 1}. We will also denote the action of the Pauli matrices and 
(T^ on the i-th qubit as Xi and Zj, respectively. 

Let A^Qj^NOT be the XOR-controUed 3-qubit gate defined on basis 1x1X2X3): 

A|orNOT( 1x1X2X3)) = 1x1X2X3) if X0R(xi,X3) = and 

A|orNOT(|xiX2X3)) = |X1N0T(X2)X3) if X0R(X1,X3) = 1, 

where X0R(xi,X3) = xi + X3 mod 2, and NOT(xi) = 1 — xi. 

Set Ui-i^i^i+i = pRy^cfi), Hi = Zi^iXiZi+i if m > 5 and Hi = Xi ii m = 3, 



Vi 



e 



^\ and not; 



\ 

isin(^) 




= AioRNOTi_i,i,i+i. 

The index iG{l,---,n-|-l}as there are (n-l-l)-qubits. The qubit encoding uses 
X^ (or X'^) particles at the two ends of the trivalent tree, and {n + 1) Yi-particles 
in the middle as depicted here. We show the case where both ends are X^: 
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X — \ \ \ \ — \ \ — X 

Xi X2 X3 X4 Xn Xn+l 

where Xj = or X[. We set Zq = Z„+2 = +1, and note that NOT^^ = / if 
Zi_iZi^i = 1 and NOT.^^ = Xi if = —1. When the proof of a following 

lemma is a direct computation, we will simply omit it. 

Lemma 3.5. 

(3.6) f/.-i,,+i = Vi ■ NOTr'+, 

Lemma 3.7. The group generated by NOT^^'^ is a subgroup of the Clifford group, 
which will be denoted as H. 

Lemma 3.8. 

(3.9) (NOTi'+y HiNOT;'^ = Hi, 

(NOTr'+yi/,+iNOTr'+ = H,H,^,, 

(N0T-'+yi7,_iN0T-'+ = H,_,H,, 



and 
(3.10) 

Definition 3.11. Define for /c < / 



i-j\>l ^ (NOTr'+)V^. 



not; 



(3.12) 
Lemma 3.13. 

(3.14) i = k-l -> 

i = k and k < I — )■ 

i = I and k < I — )■ 
i = l + l 

Otherwise — )■ 



Sk,i = Yl Hj. 

k<j<l 



[NOTr'+ 
('NOTr'+ 
('NOTr'+ 
('NOTr'+ 
('NOTr'+ 



5fc ,NOTr'+ = Sk-i 



Sk 



Sk 



:NOTr'+ = Sk+i,i, 
:NOTr'+ = Sk,i-i 
:NOTr'+ = Sk,i+u 
:NOTr'+ = Sk,. 



Theorem 2. For each n > 2 the image Pr^^ (^n) of the n-strand braid group !B„ 
is a finite group G H, where G = Z2 x and H is the subgroup in Lemma 
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Proof. The image group pi?^^ is a subgroup of the group generated by the 
generators NOT^'"^ and Vi. Note that 

(3.15) [V^,VJ]=0 
for all i,j and similarly 

(3.16) [Hi,Hj]=0 

for all Consider a word W in this (possibly larger) group consisting of a 
product of generators NOT"'"'' and Vi. Our strategy is to commute the Vi to the 
right until the word is brought into a form of a product of NOT"'"'' generators 
followed by a product of V generators (all V generators appear on the right). 

Note that NOT"'"*" is in the Clifford group. So, it conjugates ZjXj+iZj_,_2 to 
another product of Pauli gates. The matrix NOT^'"*" commutes with Hi, and 
Note that NOT"'"*" trivially commutes with Hj for |i — j| > 2. By Lemma 
13.81 NOT"'"*" conjugates a product of Hj to some other product of Hk. It follows 
that NOT"'"^ conjugates an exponential of a product of Hj to an exponential of 
a product of Hk. Therefore, any word W can be written as a product of NOT"'"'' 
followed by an exponential of a sum of terms, each term being of the form il^ 
multiplied by a product of H/s, where / is some integer. 

The group generated by the operators Vi is an Abelian group, which we call 
G. Since V^^™" = 1, we can write an arbitrary element of the group as Hj^-^'') 
where the kj G {0, 2m — 1}, so the group is a subgroup of Z2^. However, since 
Vl^ = —1, there are only 2m" distinct group elements which can be written as 
(1) Ylj Vj \ where the kj G {0, ...,m — 1}. This group is in fact Z2 x Zj!^, and the 
generators of the group can be taken to be —Vi and —1. The group generated 
by the operators NOT^'"*" is a subgroup of the Clifford group; call this group H. 
Then, because conjugation by NOT^'"*" defines an automorphism of G, the group 
generated by Vi and NOT^^'^ is the semi-direct product G y\ H. 

□ 

The finiteness of braid group images for the representations py- also can be 
deduced from |NRj . But the result in |NRj does not give information on what the 
group is. However, even our approach here does not give the complete information 
on the group. 

4. Generalized Knill-Gottesmann Theorems 

In this section, classical simulation of quantum circuits always refers to efficient 
classical simulation. 
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4.1. Classical simulation of braid group representations. Efficient classical 
simulation of Clifford circuits is provided by the Knill-Gottesman theorem (see e.g. 
[NC] ). In the context of anyonic quantum computation, it implies that braiding 
quantum circuits based on the Ising anyon a can be efficiently simulated classically. 

Fix a qudit C"^,m > 3, the state space for n-qudits is (C"^)®". A quantum 
circuit model consists of a few fixed unitary matrices {gi}, called a gate set, and 
all quantum circuits on n-qudits based on this gate set for all possible n > 1. A 
quantum circuit on n-qudit is a composition of finitely many gates in the sense 
that gates should always be extended by tensoring an appropriate identity matrix 
if necessary. 

Given a i?-matrix : C™ ® C™ — ^ C™ ® C™, by assigning 

/^(^-^) ^Rx(S) J^(""^-i) 

to the braid generator (Tj, z = 1, ■ ■ ■ , n — 1, we obtain a braid group representation 
of ^n- Naturally, we may regard this as a quantum circuit model with the 
gate set {pr^ (^i^)}' — 1)2,- ■■ ,n — 1. Then for each braid a E S„, pn^i^a) is 
a braiding quantum circuit on n-qudits (C™)*^" for all n > 2. Note there are no 
1-qudit gates. 

Conjecture 4.1. If a unitary (generalized) R-matrix Rx is of finite order, then 
braiding quantum circuits can always he simulated classically. 

In this section, we will prove that this indeed is the case for the localizing qudit 
representation p^^ , and also the generalized localizing qubit representation p/j^^ . 

4.1.1. Ising anyon a. The Ising anyon a can be localized by the following i?-matrix 
[FRWj . which consists of the Bell states. 





/ 1 








1\ 


1 





1 


-1 





7! 





1 


1 















V 



As a corollary of the Knill-Gottesman theorem, any braiding quantum circuits 
from this i?-matrix can be simulated classically. 

4.1.2. Localizing the Gaussian representation. We want to find a unitary m? x m? 
matrix U so that Ui := Id®^-^ ®U ® Id®'""-'-^ satisfy (13.1113.31) . For this let {ei}T=i 
denote the standard basis for C" and define 

(4.2) U (ej ®ej) = uj'~^ei+i ® e^+i 

where the indices on e^ are to be taken modulo m. It is straightforward to check 
that t/™ = / and U'^ = U^^ . It is clear that Ui and Uj commute if |i — j| > 1. It 
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remains to check (13. 2p . For this it is enough to consider i = 1: 

UiU2{ei (g) (g) efc) = a;*~''"^ei+i (g) e^+s ® e^+i = 

^_2^i-fc+lg^^^ (g) ej-+2 ® Bfc+i = W"^f/2f/l(ei (g Bj (g Gfc) 

so lij — t/j does give a representation of ES{u, n — 1). A standard trace argument 
shows that this representation is faithfuL Thus defining 

^ m— 1 
/m < ' 



m 

j=0 

gives an i^^^-matrix locahzing the Gaussian representation, therefore the braid 
group representation px^,n |RWt IRoj . The resulting braid group representation is 
the qudit representation px^, denoted as Pr^^- 

4.1.3. Qudit representation pR^^ . For the qudit space C™, we denote its standard 
basis as {|j),0 < j < m — 1}. On the i-th qudit C", we define the "shift" 
and "clock" operators Xj and Zi. The "shift operator" Xj is the permutation 
matrix Xj|j) = \j + 1), and the "clock operator" Zi is the diagonal matrix Zi\j) = 
< j < m — 1. When m = 2, Xj, Zj are the Pauli matrices. 
The qudit representation is a direct sum of many sectors with multiplicities. 
Each sector has a graphical calculus using the following tree, where the two hor- 
izontal ends of the trivalent tree are labeled by two anyons Yl,Yr in {Yi,i = 
0, 1, r}, m = 2r + 1: 

Xf Xf Xg Xg Xf Xf 

Yl—\ \ \ L \ Ly^ 

Xi X2 ^n-l 

Theorem 3. When m is a prime p, then all polynomial length braiding circuits 
of anyons can be efficiently simulated classically in the sense of the stabilizer 
formalism. 

Proof. There are (r + 1)^ different sectors. Since braiding X^ anyons will not mix 
sectors, so they can be simulated simultaneously or individually. The representa- 
tion space of px^,n is Hom(yi(g (X^)®", Yr) with a basis consisting of eigenstates of 
products of Mj's. They are labeled trees in Fig. 14.1.31 By equation (13. 4p . braiding 
Xe anyons transforms products of Mj's into products of Mj's. As a result, the evo- 
lution of a state in B.om(Y]^ (g (XJ*^", Yr) can be efficiently simulated classically 
by following the evolution of these operators. 

The extra-special p-group 2k another generating set [HNW] : 

XjXj = XjXi , ZiZj = ZjZi 
XiZj = z ZjXi 
(4.3) XiZ = zXi , ZiZ = zZi 
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These two presentations of -E^ 2fc are related by M2i-i = -^i, M2i = ^i-^i+i ior i ^ k 
and it2fe = Zh- We will include E^^ ^k-i inside -Ep 2fc- 

The extra-special p-group -Ep 2fc introduces redundant states for X^-anyons, which 
will be removed by the stabilizer formalism. Braidings commute with roughly half 
of the generators of E^^k- Note that Ui, ...Un-i, Ui, ...Un^i, XiZi, XnZ^, z, where 
Ui = XiXi^iZiZj_^-^_,Ui = XiXi^iZjZi^i is another generating set of Ep^2k- The 
generators Ui,XiZi, and XnZl all commute with the Ui's and, therefore, with 
braiding. 

It suffices to show that the two-qubit gates Rx^ conjugates each of the gen- 
erators Xi,X2,Zi, Z2 to a product of these generators, up to a phase. A direct 
computation shows: 



(4.4) U = XiX2ZiZ2\ 



then U commutes with the operators X1X2, Z\Z\^ X-^Zi, and X2ZI as may be 
checked. So, all these 4 operators are mapped to a product of the generators. 
Note that these 4 operators are not all independent: the product of the first three 
operators is equal to the fourth operator, up to a phase. 

Note that the operators Zi, X1X2, ZiZl, X^Zi generate the image group. Hence, 
it suffices to check that Zi is mapped to a product of generators by Rx^, up to 
phase. Let I — Note that 



(4.5) 



WZi = uj-^ZiW. 
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Then we have: 

(4.6) RxZ,R\^ 

p-i 

i=o 

p-i 

i=o 
i=o 

i=o 

p-i 
i=o 

Therefore, the evolution of Zi can be efficiently simulated classically and, con- 
sequently, so can the evolution of any state stabilized by products of Xi and Zj 
operators. □ 

For the standard quantum circuit model, measurements in the middle of com- 
putation can be all postponed to the end (see e.g. |NC] ). But this is not the case 
for the anyonic computational model. So we may ask if we can also simulate some 
measurements in the middle of a computation. Measurements in a basis of prod- 
ucts of Xi and Zj operators can be simulated classically as in the Clifford circuit 
case. Thus, we conclude that we can efficiently simulate classically any polynomial 
length quantum operation that consists of creating pairs of X^ anyons out of the 
vacuum, braiding them, and then measuring them in a basis of products of Xi and 
Zj operators. 

More difficult measurements are projections of a pair of X^ anyons onto a definite 
charge. While projecting onto the trivial charge can be simulated classically, we 
do not know if this is true for nontrivial charges. This question seems to be open 
even in the Ising case. 
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When m is not a prime, then the equivalence of the Gaussian representation 
and the braid group representation px^ is conjectured to be true, but not known. 
Note that the exact same proof gives a classical simulation of the Gaussian repre- 
sentation. 

Instead of u"^ = 1, if we set wf™ = 1, then the Gaussian representation and 
its localization are defined similarly. The same proof above will give a classical 
simulation of the localizing braid representation. 

4.2. Qubit representation p/j^^ . If Ry : V®^ — )■ F®^, V = C™ is a solution to a 
(2, 3, l)-gYBE, for the braid generator CTj G 23„, by setting 

we again have a braid group representation |KWj . But in this case, !B„ acts on 
the vector space The above discussion about braiding quantum circuits 

applies to the generalized localizations as well. 
We will use the same notation as in Section 13. 3[ 

Theorem 4. All polynomial length braiding circuits from p^^^ of S„ after a gen- 
eralized localization can he efficiently simulated classically in the sense that we can 
efficiently classically represent and multiply elements ofGxiH. 

Proof. As explained above, we can simulate each sector independently. In Section 
13.31 the image group is given as a subgroup of a semi-direct product of two groups 
G xi H. G is an abelian group with generators e* ™ for each pair of integers k, I 
with l<A;</<n + l. Thus, the group G has O(n^) generators and is finite 
order so we can write each group element by writing a sequence of integers in 
the range {0, 2m — 1}. 

The other group H is generated by the NOT^^'^ operators. This is a subgroup 
of the Clifford group. For a unitary U in the Clifford group, let us compute 
UXiW and UZiU'^ for all i. If one knows these conjugations of Pauli operators, 
then one can compute UOW for any operator O since any O can be written as 
sums of products of the Xi and Zi, where Xi and Zi are Pauli 's. Since a unitary in 
the Clifford group conjugates Pauli's to products of Pauli's, we can store UXiW 
and UZiW for all i with polynomial resources. 

So, this describes how we will represent elements of the two groups G and H. 
Note that we only are storing U up to an overall phase. The classical simulation 
is done as follows. After some number of steps, we have a unitary U stored as 
AC, A & G,C & H, where we store A by storing the integers and we store C by 
its action on the Pauli's. We want to left-multiply f/ by a generator which is either 
a generator in the Abelian group G or an element of the group H. Note that the 
braid group representation image is generated by a product of an abelian group el- 
ement A' and a G' operator, but we describe the left-multiplication by the abelian 
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elements A' and the C operators separately. Left-multiplying by an abelian gener- 
ator A' is easy: simply multiply A by the abelian generator. Left multiplying by a 
C is done in two steps by considering the identity C'AC = {C'AC'^)C'C = A'C'C. 
First, we commute the C through the element A. Using the commutation rela- 
tions we have given, we can do this multiplication in polynomial time: each S^^i 
gets conjugated by the C to some other Sk^i, so commuting the C through A 
permutes the different integers describing A. Next, we left-multiply C by the C; 
we can do this by taking CXjC^ and conjugating by the C which can be done in 
polynomial time, and similarly for CZiCK 

□ 

Since any operator Ski is conjugated by this group to another Ski, we are able 
to classically simulate braiding combined with measurement of the operators Ski- 
Physically, the measurement of Sj^j+i corresponds to whether or not the i-th and 
{i + l)-th anyons fuse to a Yi anyon or to 1 or Z. However, this measurement 
does not distinguish the fusion outcomes 1 or Z from each other. Ultimately, this 
is related to the #P-hardness of evaluating the related link invariant in the next 
section. 

5. #P-HARDNESS OF EVALUATING LiNK INVARIANTS 

In this section, we relate the evaluation of the link invariant Iyi{L) for certain 
links to the computation of the Ising model partition function of some graphs G. 
By our construction, the maximum cuts of the graphs G correspond to ground 
states of the Ising model. The evaluation of the link invariant Jy^ (L) would im- 
ply the counting of the maximum cuts of the graphs G, a well-known #P-hard 
problem. 

Suppose J{L]x,y) is a polynomial invariant of oriented links and po = (xo,?/o) 
is an algebraic point that all evaluations J(L;po) exist. Then J(L;po) is a numeri- 
cal invariant of oriented links, which sometimes can be identified with a partition 
function of a (2 -|- 1)-TQFT. We are interested in the computational complexity 
of such evaluations po) for all oriented hnks. To be precise, we need to spec- 
ify an input encoding for links. There are several equivalent encodings such as 
combinatorial data for link diagrams or words in braid generators. 

There are two 2-variable generalizations of the Jones polynomial: the HOMELY 
polynomial P{L; I, m), and the Kauffman polynomial F{L; a, z). The Jones poly- 
nomial V{L] t) is the specialization of the Kauffman polynomial F{L\ t~^, — (t~4 + 
t*)). Eor these two polynomials, the complexity of evaluating each algebraic point 
has been shown to be either ^P-hard or polynomially computable classically |We] . 
When the computation is polynomial classically, the evaluation has an interpreta- 
tion using classical topological invariants. We will refer to the classically polyno- 
mial time computable points as the classical points. 
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5.1. Evaluating Ix^{L). For the simple object in S0{m)2, the resulting link 
invariant Jx, [L) has a classical interpretation. They are not specializations of the 
HOMFLY or the Kauffman polynomials when m > 7, and satisfy a (^^Y^)-term 
skein relation. When m = 3, IxX^) is the Jones polynomial at a 6th root of unity 
and when m = 5, an evaluation of the Kauffman polynomial at a 5th root of unity 
[Jo89] . 

Proposition 5.1. The exact evaluation of Ix^L) is polynomial time classically 
when m is a prime p. 

In |Jo89t IGJoj . Ix^L) is expressed as a sum using the symmetrized Seifert 
surface 5 of a braid closure: 



^ veHi{S;¥p) 



, ,<v,v> 



where <, > is the Seifert form. The norm of this sum is (i/p)^, where r is the rank 
of the first mod-p-homology of the 2-fold branched cover of the 3-sphere along 
L. The phase is given by a Legendre symbol. 

5.2. Evaluating Jy^ (L) . For the simple object Yi of JKr=Q and its generalizations, 
the resulting link invariant normalized to oriented links is essentially the evaluation 
of the Kauffman polynomial at (a, = (— i e~^, 2sm(^)). Precisely, 

/y,(L) = i(-irW-iF(L;-ie-^,2.m(-)) 
/ m 

for a c(L)-component oriented link L |Hong| . When m = 3, Iyi{L) is the evaluation 

of the colored Jones polynomial at a 6th root of unity. 

Theorem 5. The exact evaluation of ly^L) is ^P-hard. 

Proof. For the Kauffman polynomial F{L; a, z) such that a 7^ 0, 2; 7^ 0, the classical 
points are |Wej : 

(1) a = ± i; or 

(2) (a, z) = {—q^^, q + q^^), where q^^ = 1 or g^^=l, but g 7^ ± i; or 

(3) (a, z) = (g^^, q + q~^), where g*^ = 1 or g^^=l, but g 7^ ± i; or 

(4) (a, z) = (— g^\ g + g~^), where g^^ = 1, but g 7^ ± i; or 

(5) (a, z) = ±(1, g + q'^), where g^ = 1. 

Noticing that (a, z) = (— ie~^, 2sm(^)) are not in the list, we obtain the 
desired result. 

□ 

This proof of the #P-hardness of evaluating Iyi{L) is very indirect. To have a 
better understanding of the #P-hardness, we provide two refinements. 

Theorem 6. (1) Any sufficiently accurate approximation of the link invariants 
IyAL) ts i^P-hard. 
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(2) The sign of ly^ (L) is jj^P-hard to approximate. 

The specialization of the Kauffman polynomial F{L\ a, z) to z = — (a + a~^) has 
a state sum due to R. Lickorish and K. Millet |Lic] : 

E{L) = 2(-l)^(^)-iF(L) = ^a-^<^'^-^\ 

S(ZL 

where the summation is over the 2'^'^^^ sublinks of L including the empty link, 
a = —i exp(— ivr/m), and {S, L — S) is the total linking number of the two sublinks 
S and L — S. This state sum of Iyi{L) allows us to translate Ising model partition 
functions to Jy^ (L) for certain links. 

5.3. #P-hardness of approximation. We first define a partition function for 
an Ising spin system, a quantity frequently used in physics and closely related to 
the Tutte polynomial. The restriction to integer entries for the matrix below is 
because those will be the only cases we need. 

Definition 5.2. Let J = {Jij)i<i,j<N be a symmetric N x N matrix with integer 
entries and with diagonal entries of J equal to 0. Let y be any real number. Define 
the Ising partition function to be 

(5.3) Z{J,y) = Yl y^^<^''''^-''^, 

ae{-l,l}'^ 

where the sum is over A^-dimensional vectors a with entries chosen from the set 
{ — 1, 1} and where S^j.^^. is the Kronecker symbol. 

The quantities y^i<J '^'•'^^ are called the Boltzmann weights. 

Often in physics one restricts to y > so that the Boltzmann weights are 
positive, but we will not make this restriction here. Also, often in physics one 
considers a closely related quantity: 

(5.4) ^ exp(-/3^ JijCTiCTj). 

cre{-l,l}^ i<3 

Setting 



(5.5) 



P = -\Hy). 
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the quantity of Eq. (15 .4^ is equal to 

J2 exp(-/3 5^J,,(25.,^,.^,-l)) 
Z{J,y)exp{/3^Jij) 

i<j 

so that Eq. (15. 4p is the same as Eq. (15. 3p up to a trivial multiplication by exp(/3 X]j<j "^ij)- 
We prefer to use Eq. (15.31) because it will be simpler later when dealing with neg- 
ative y. 

Lemma [5?7I below shows that for certain J, y, we can define a link whose link in- 
variant is equal to the partition function Z{J, y), up to multiplication by a quantity 
that can be computed easily. The allowed y take certain discrete values depending 
upon m. The allowed y obey \y\ < 1, and we will choose y such that \y\ < 1 to ob- 
tain a nontrivial Z{J, y). Importantly, the size of the link will only be polynomially 
large. 

We then give two different proofs of #P-hardness, each leading to a slightly 
different improvement of the T^P-hardness. In the first proof, we choose y and 
J such that all the Boltzmann weights are positive; then, by taking the coupling 
constants in the Ising system large, a sufficiently accurate evaluation of the Ising 
partition function allows one to count the number of maximum cuts a given graph 
has, a known ^P-hard problem. The second proof is based on considering a 
different regime using negative Boltzmann weights; in this known result 

|GJj shows that evaluating the sign of Z{J,y) is #P-hard. 

The sign of the Boltzmann weights depends upon the sign of y and on the 
parity of the entries of J. For every odd m > 3, there exists a choice of d, with 
1 < d < m — 1, such that y in Eq. (15. 8 p below is negative and greater than —1; 
this regime will give a complex /3. In this case, the Boltzmann weights may be 
negative if J has odd entries and this case is the basis of the second proof. To 
obtain positive Boltzmann weights, we can either restrict to matrices J with even 
entries, or we can consider odd m > 3, for which it is possible to pick d such that 
y is positive and less than 1. To keep a unified treatment for all m, we follow the 
choice of using matrices with even entries in the first proof. 

The positive Boltzmann weight approach has the advantage that the quantity 
that we show is T^P-hard to approximate has a simpler experimental realization 
than the second regime does: the absolute value squared of this quantity is the 
probability that preparing pairs of Y particles that fuse to the vacuum, then braid- 
ing them in a way determined by the link, and then fusing them together in pairs 
will have all pairs fuse to the vacuum, while the sign is not in principle measurable 
without doing an experiment that interferes different particle trajectories. Also, 



(5.6) ^ exp(-/3^ Jijo-jO-j) = 

o-e{-i,i}^ i<i 
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the result below that it is #P-hard to approximate this quantity to a multiplica- 
tive accuracy better than a given amount (this accuracy is exp(— poly(A^)) can be 
straightforwardly translated into a result that it is ^P-hard to approximate to an 
additive accuracy better than a given amount since we can give a lower bound on 
the quantity itself, being a sum of positive terms. 

The negative Boltzmann weight approach of lemma 15.331 below has the ad- 
vantage that since we show that it is ^P-hard to approximate the sign of the 
link invariant, it shows that approximating the link invariant to any (positive) 
multiplicative accuracy is hard, improving greatly on the exponentially small mul- 
tiplicative accuracy required in the first regime. 

Lemma 5.7. Let J = (Jij) be a symmetric N x N matrix with Jij G Z, Jjj = 0, 
and d be any integer. Set 

(5.8) y= ^ , 

•where 

(5.9) a = — i exp(— iyr/m). 

Define P{J) to be the sum of the positive entries of J, and A{J) = J2i<ij<N 
Note that P{J) and A{J) both are even since J is symmetric. 

Then, there is a link L containing 0{N'^) components such that the quantity 
E{L) = 2{-l)<^^'^F{L) defined by Lickorish and Millett [Oc] obeys 

(5.10) E{L) = Z{J, y)a-^^^(^) (y^) " ( v/2(a-4^ + a+^^)) 

where c{L) is the number of components of the link and where we choose the sign 
of the square-roots such that y^-\/2(a~^^M-a+^ = a~^'^ + a~^^^. 



The number of crossings in the link is at most polynomial in N + dA(.J). 
Proof. Lickorish and Millett show that E{L) can be written as a sum 

(5.11) E(L) = ^a-^<^'^-^\ 

ScL 

If a link L has c(L) components, then a sublink S can be specified by specifying, for 
each component, whether that component is in 5* or not. We do this by defining, 
for each sublink S, a vector s with entries Si for i = 1, c(L), such that Si = +1 
if the z-th component is in S and Si = —1 otherwise. Let {i,j) be the link number 
between the i-th component and the j-th component. The invariant {S, L — S) is 
then equal to the sum of {i,j) over pairs i G S and j E L — S. 
So, 

(5.12) - 4(^, L - S) = - 5^(1 + s,)il - s,){t,j) = -2 ^(1 - SiS,){t,j) 
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So, Eq. (15. lip is equal to 

(5.13) E{L)= J2 a^^^<^'^'''^''^^a-^^'<^^''^l 

se{-i,i}':(^) 

We now show how to construct the desired hnk L. We pick 

(5.14) ciL) = N + J2\J^j\- 

i<j 

So, we can define a one-to-one mapping from the set of triples of integers, ^, j, n 
with 1 < i < j < N and 1 < n < \Jij\ to the set of integers from + 1 to c{L). 
We describe this mapping by a function F{i,j,n). 

Pick = if both i and j are less than or equal to or if both i and j are 
greater than or equal to A^. Consider each triple i,j,n as above. Then pick 

(5.15) {t,F{t,j,n)) = d 
for all n. Pick 

(5.16) (j, F{i,j, n)) = d- sgn{Jij) 

for all n, where sgn(a;) is the sign function: +1 for a; > and —1 for a; < 0. 
Note that by definition {i, F{i, j,n)) = {F{i, j,n),i) and similarly {j, F{i, j,n)) = 
{F{i,j,n),j). Let all other link numbers {k,l) = 0. 

This completes the description of the link; we now show that the link invariant 
E{L) is the desired result. We do this by defining 

Si=ai,l<i<N 

(5.17) X(L,a)= g 

se{-i,i}=(^) 



The summation notation means the sum over all c(L)-dimensional vectors s with 
with entries +1 or —1, such that the first entries of s are equal to those of a 
for some A^- dimensional vector a. Then, 

(5.18) E{L)= X{L,a). 

ae{-i,i}^ 

To compute X(L, o"), note that this is equal to 

\Jij\ 

(5.19) JJ^ ( (x^'^i'^F(i,j,,^){i,F{hhrl))+2sjSp^^^j^^){j,F{ij,n))'y^ 
l<i<j<n n=l SF(,,,,„)e{-l,l} 

l<i<j<n 
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where 



\Jro\ 



(5.20) n( E 



W{i,j,Jij) = 

ij,„)(i.f'(i.i.")>) 



2siS 



i^F(i,j,n) (i,F(i,j,n))+2s:s 



n=l s 



a' 



■5F(i,j,n)G{-l.l} 



A direct computation gives in the case Sj = sgn{Jij)sj that W{si, Sj, Jijj 
{a~^'^ + a+'^°')l'^^jl and in the case that Si = —sgn{Jij)sj that W{si, Sj, Jij) = 2^-^'^^. 
Let y be as defined in Eq. (15.81) and let 

(5.21) z = 2(a-^'^ + a+^'^). 
Then 

(5.22) W{s,, s„ J,,) = v/^-'-'^'^^V^'^-', 

where the ambiguity in the sign of the square-root is resolved by choosing ^/y^/z = 
So, 

(5.23) X{L, a) = y^^<^ ^-'^^-^ (v^) ^"a-2'^^(-^)v^^^<^ I-"-'. 

Comparing to Eq. f l5.10p completes the proof. □ 

Definition 5.24. Given a graph G, with a set of vertices V , define a cut of G 
to be a partitioning of the vertices of the graph into two subsets, 5* and V — S. 
Define the size of a cut 5, \^ — to be the number of edges in the graph connecting 
vertices in S to vertices inV — S, and denote it as size(G,S). Define M{G) to 
be the maximum over all cuts of graph G. Define N{G) to be the number of cuts 
with size equal to M(G), so that 

(5.25) N{G) = \{S<Z r|size(G, S) = M(G)}|. 

Note that according to this defintion we consider S", V — 5 to be a distinct cut 
from V — S,S, so that a cut is defined by an ordered pair of sets. 

Lemma 5.26. Consider any N x N matrix J with all diagonal entries equal to 
zero and with all other entries equal to either or K for some even integer K . 

Define a graph G to have N vertices, labelled by integers between 1 and N, and 
to have edges connecting vertices i,j if and only if Jij ^ 0. 

Fix d such that Eq. I{5.8\) gives \y\ < 1. 

Then, the quantity Z{J, y) defined in Eq. h5. 3\) obeys 

(5.27) iV(G)|i/|^(l^l-^^(^)) < Z{J,y) < iV(G)||/|^(l^l-^^(^» + \y\K{i+\E\-MiG))^N ^ 
where \E\ is the cardinality of the edge-set of the graph so that K\E\ = Jij- 
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Proof. There is an obvious one-to-one mapping between choices of the vector a in 
Eq. f l5.3p and cuts: for a given choice of the vector a, define a set S to contain all 
vertices i such that o"j = +1 and define — to contain all other vertices. Then, 
Eq. (15. 3p gives Z{J,y) as a sum over cuts. For a maximum cut, the Boltzmann 
weight equals \y\^^\^\~^^^'^')\ while for all other cuts the Boltzmann weight is at 
least a factor of \y\ smaller. Since there are only 2^ cuts and since every cut 
contributes a positive term in the sum, the desired result follows. □ 

Lemma 5.28. Fix d such that Eq. \5.3^) gives |y| < 1. 

We consider symmetric matrices J such that all diagonal entries are equal to 
zero and all other entries are equal to either Q or K for K chosen to he the smallest 
even integer greater than or equal to 

(5.29) (iV + l)ln(2)/ln(|t/|). 

Consider the problem of approximating Z{J, y) to a multiplicative accuracy within 
the range [1 — 2^^~^, 1 + 2^^^^]. That is, we want to compute a number Z such 
that it IS guaranteed that (1 - 2~^-^)Z < Z{J,y) < (1 + 2-^'^)Z. 

This approximation problem is jj^P-hard. 

Proof. Consider an arbitrary graph G with vertices. Let the matrix J he K 
times the adjacency matrix of G. We will show that being able to evaluate Z{J, y) 
to the given accuracy allows one to compute M{G) and N{G). 

Note that N{G) < 2^. Also note that for the given K, \y\Kii+\E\-M(G))2N < 
(l/2)|?/|^(l^l-^(^)). Therefore, all graphs G with a given M{G) have a Z{J, y) that 
is at least twice as large as that of any graph G' with M{G') < M{G). So, given 
an approximation Z with the given accuracy, it is possible to determine M{G) 
exactly since M{G) is an integer. 

Next we show that we can determine N{G). For the given K, Eq. (I5.27P gives 

Ar(G)|y|^(l^l-*-^(^)) < ZiJ,y) < iV(G)|y|^(l^l-"-^(^» + \y\m+\E\-AHG))2N 

< (Ar(G') + l/2)||/|^(l^l-^'^(^)). 



Equivalent ly. 



Z{J,y) 



(5-30) NiG) < ^y^^^;,;_^^a)^ < NiG) + 1/2. 

Let 

(5.31) N 



K{\E\-M{G))- 



\y 



So, (1 - 2-^"3)iV < N{G) + 1/2 so 

(5.32) (1 - 2-^-'^)N - 1/2 < N{G) < (1 + 2-^~^)N. 
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Since A^(G) < 2^, it follows that iV < (2^ + l/2)(l+2-^-3) < 2^+1 So, < 
1/4. So, since N{G) is an integer, the above equation determines N{G). 

So, given this approximation, it is possible to determine both M(G) and N{G). 
That is, it is possible to determine both the MAX-CUT of a graph (an A^P-hard 
problem) and to count the number of such maximal cuts, a ^P-hard problem 
|Sharp| . □ 

Lemma 5.33. Fix d such that Eq. / 1 5. <g)) gives —1 < y < 0. Then, the quantity 
Z{J,y) is real but it is ^P-hard to compute its sign. 

Proof. Choose the entries of J to be either or 1. Then, a computation of the 
sign of Z(J,y) gives a computation of the sign of the Tutte polynomial Tci^x^y) of 
a graph G whose adjacency matrix is J with y as given here and with x such that 
(x — l){y — 1) = 2. However, computing this sign is #P-hard |GJ] . □ 

Our improvement for #P-hardness is to relate the link invariant Iyi{L) to the 
Ising model partition function. For certain links, the link invariants allow one to 
approximate the Ising partition function. We have seen that approximating the 
Ising model partition function has at least two different regimes, one regime where 
it is ^P-hard and one regime where it is A^P-hard, depending upon whether we 
compute both N{G) and M{G) or just M{G). Approximating Iyi{L) is more 
general than just approximating the Ising model partition function at real tem- 
perature because the corresponding temperature could be complex, i.e., by picking 
some very specific links, we constructed a problem with a real temperature, but 
for arbitrary links we have to deal with a complex temperature. However, if we 
assume that all we are doing is approximating the Ising model partition function 
when we approximate the link invariant, the problem goes from A^P to t^P as we 
approximate it more accurately, while link invariants of computationally universal 
theories like SU{2)^ instead go from BQP to #P, which is quite different. Since 
BQP is believed to be incomparable to A^P, there is no approximation regime in 
which such other link invariants would be A^P. It is not clear if this A^P-hard 
regime will survive complex temperature or not, though, but it would be a very 
interesting phenomenon. 

In |KRj . the quantum double of finite groups are studied. Some theories could 
also have finite braid group images but #P-hard link invariants. It is likely that our 
theory follows the same pattern of their analysis, where different approximations 
to the link invariant are in P, or are S'PP-hard, or are #P-hard, depending upon 
the accuracy of the approximation. It would be interesting to see if quantum 
doubles can be efficiently simulated classically. 
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